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Abstract: We define a new multivariate time series model by generalizing the ARMAX process in a multivariate 
way. We give conditions on stationarity and analyze local dependence and domains of attraction. As a conse- 
quence of the obtained result, we derive a new method of construction of multivariate extreme value copulas. 
We characterize the extremal dependence by computing the multivariate extremal index and bivariate upper tail 
(*~*) dependence coefficients. An estimation procedure for the multivariate extremal index shall be presented. We also 

£Nj , address the marginal estimation and propose a new estimator for the ARMAX autoregressive parameter. 
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Ov 

E— i ! 1 Introduction 

Stationary time series presenting sudden large peaks are usually well modeled by heavy tailed noise 
ARMA. However, models with practical application but simpler treatment have been studied in literature 
as an alternative. Davis and Resnick ([5], 1989) proposed the MARMA process which is analogous to the 
ARMA by just replacing summation by the maximum operation: 



X, = faXi-! V ... V 0pX;_ p V Yi V OiYi-! V ... V 6 q Yi- q , 

>' 

" where < <pi,0j < 1, 1 < i < p, 1 < j < q and the innovations Y n , n > 1, are independent with unit 

OO ■ Frechet distribution. A first order MARMA type version, i.e., 

00 



(N 



X i = cX i _ 1 VY i , 0<c< 1, (1) 



was analyzed in Alpuim ([I], 1989) by considering innovations Y n , n > 1, independent and equally 
£\| . distributed, not necessarily unit Frechet. The model in ([T]), sometimes denoted in literature as ARMAX, 

corresponds to the case a = in the Haslett ([H], 1979) model 



X i = 0X i - 1 V(a0X i - 1 +Y i ), 

X. 

0<p<l,0<a<l, used to describe a solar thermal energy storage system and later developed in, 
e.g., Daley and Haslett ([2] 1982) and Greenwood and Hooghiemstra ([121, 1988). Exponent versions of 
ARMAX, namely pARMAX and pRARMAX, were used in the modeling of financial series (Ferreira and 
Canto e Castro, [TT] 2010). Further applications of ARMAX processes and their generalizations can be 
seen, for instance, in Lebedev ([T7], 2008) and references therein. Here we shall consider a multivariate 
formulation of ARMAX and extend some of the results in Alpuim (pQ, 1989). More precisely, we will 
analyze conditions on stationarity (Section [5]), local dependence conditions (Section [3J and domains of 
attraction (Section 0]). The relation between the max-attractors of the process and the innovations allow 
us to evidence a new construction method of copulas for multivariate extreme value distributions (MEV) . 
In computing the multivariate extremal index, we shall find that it is possible to have clustering in all 
or only in some of the marginals, according to their domains of attraction. An estimation procedure 
for the multivariate extremal index will be also stated (Section 0|. In Section [S] we will derive the lag-r 
tail dependence coefficient (TDC) (Sibyua [53], 1960; Joe [T3], 1997) and the lag-r tail independence 
coefficient of Ledford and Tawn ([HHH], 1996/97) and we find different types of tail dependence. Some 
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notes on the marginal parameters estimation shall be given at the end (Section 12]). In particular, we 
present a new estimator for the ARMAX parameter c which is strongly consistent and asymptotically 
normal. 



2 The multivariate model 

Let {X„ = (X nt i, X n ^)} n >i be a d-variate sequence, such that 

X n . 3 = CjX n _ x<j V Y nJ , n > 1, j = 1, d, < Cj < 1, (2) 

where X = (X ,i, X 0td ), Y„ = (Y nA , ...,Y n , d ), n > 1, are independent, X — F and Y„ ~ G. 
{X„}„>i thus corresponds to a d-variate formulation of an ARMAX process given in ([1]). 

It is an immediate consequence of relation ([2]) that each marginal {A„j}„>i, j G D = {1, ...ci}, of the 
sequence {X„}„>i can be written in the form 

n 

X n>j = c]X 0j V \/ c]-%j, n>\. (3) 

i=l 

If F n denotes the distribution of X„ = (X n i, ...,X nt d), we also have 

F n (x 1} ...,x d ) = F„_i — , — )G(x 1 ,...,Xd) (4) 

V Ci C<j / 

and 

n 

F n (x lt ...,x d ) = F (j|, 5jp II G (^I=T' -' ^ 

If there exists lim^-^oo F n (x\ 7 then, based on we have 

— , — )G(xi, ...,Xd). 
ci c d / 

and, by ©, 



(• r >) 



(6) 



n oo 

F( X1 , *<) = Inn n Gft, .... ^) = II ^ ^ ■ ( 7 ) 

i— 1 L a i—l ± a 



Proposition 2.1 {X n }„>x is a stationary sequence with common non-degenerate distribution if and only 
if there exists (x±, ...,Xd) € such that 



0<^-logG(| J ...,^)< O c. (8) 
In this case, the common distribution F o/{X n }„>i satisfies 

Proof By 0, F is non degenerate if and only if there exists x = {x\, Xd) such that < F(x\, Xd) < 
1, i.e., such that 

0<fiG(£.™. *?)<!• 

The assertion in ((SJ) is straightforward by taking logarithms and, if it holds for some x G M d , then 
g(&,...,^iO ->■ 1, and thus x G El. □ 
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As a consequence of (|HJ|, if any of the marginals Gj, j € D, of G has support with non positive right 
end-point then the corresponding marginal F n j has degenerate limiting distribution and, therefore, F is 
a d-dimensional degenerate distribution. Observe that © is satisfied by every multivariate distribution 
with positive dependence and marginals Gj satisfying < ~^°S^j( x / c )) < 00 f° r some x > 0. 

This latter condition is satisfied, for instance, by the Generalized Pareto distribution (Alpuim, [1] 1989). 

For each j £ D, suppose that Fj belongs to the max-domain of attraction of Hj, in short Fj £ T>(Hj), 
i.e., there exists constants {a n j > 0} n >i and {b n j} n >i, such that 

n(l - F(a nJ x + b nJ )) -> - log Hj (x) , 

n— »oo 

where may be a Gumbel, a Weibull or a Frechet distribution, respectively, A(x) — cxp(— e~ x ), 
^ a Ax) = e~(~ x ~) 3 , x < 0, and $ a . (x) = e~ x 3 , x > 0, for some <x, > 0. Therefore, a sequence of 
normalized levels {uj 3 j} n >i for {X„j} n >x, i.e., such that 

n(l-^(^)) -> r,->0 

can be written as 

with Hj (x) = inf{y : F'(y) > a;} the generalized inverse of Hj. By applying the Khintchine's types 
theorem (see, e.g., Leadbetter et ai, Q~5] 1983), we arrive at the following property of the normalized 
levels for {X n j} n >i that shall be used latter: 



U n,j 



u } , with r- = <^ (9) 



In the sequel we shall denote {u„ = (u^\ , ■■■i u nd)}n>i the sequence of normalized random vectors. 

3 Asymptotic independence and local dependence of {X n } n >i 

As showed in Alpuim (pQ, 1989) for the univariate case, we shall prove that the strong- mixing condition 
also holds for the multivariate sequence, i.e., for any A £ £>(Xi, X p ) and B £ Z?(X p + s +i, X p+s+ 2 7 •■•)> 

\P(AnB)-P(A)P(B)\ <a s 



with a s —> 0, where B(-) denotes the cr-field generated by the indicated random vectors. 

s— too 

In what follows, all operations and inequalities between vectors are understood to be componentwise. 



Proposition 3.1 {X n }„>i satisfies the strong-mixing condition. 
Proof Consider A € B(Xi, ...,X p ) and B £ B(X p+s+ i, X p+s+2 , ...) and let 

C s = {Y p+ i ^ X p , Y p+s+ i ^ X p+S }. 

We have 

\P(Ar\B)-P(A)P(B)\ 

= \P{A n b n c s ) + P{A n B n c s ) - p{A)p(b n c s ) - p(A)p(b n c s )\ 
< \P(A nfinc s ) - p(A)p(b n c s )| + nsn c s ) - p(A)p(b n c s )| 
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Observe that, for the first term, 

\P(A DBDCs)- P{A)P{B n C s )\ 
= \P(C S )P(B\C S )P(A\B n C a ) - P(A)P(C S )P(B\C S )\ 

< P(C S )P(B\C S )\P(A\B n C s ) - P(B\C S )\ 

< P(C S ). 
On the other hand, since 

p+s+l-fc 

X p+S+1 = C P+ s+1 - k X k V \/ cP+ s+1 - k -' l Y k+l , k=p,...,p + s, 

we can write 

A n B n c s = A n b' n c s and BnC s = B'n c s 

where P' G $(Y p +i, Y p+ 2, •■•)■ Thus being, for the second term, we have 

\P(A DBD C s ) - P(A)P(B n C s )\ = \P(A flB'fl C s ) - P{A)P{B' n C s )\ 

< \P{A DBDCs)- P{A n B')\ + \P(A n B') - P(A)P(B n C a )| 
= P(A n B')P(C,|i4 n B') + \P(A)P(B') - P(A)P{B n C s )| 

< 2P(C S ). 

Now we just need to prove that P(C S ) — > 0. Observe that 

P(C S ) = 1 - P({J°+l{Y p+i > X p+i -i}) < f - P(US{Y P+1 > 0-%}) 

< 1 - P(Y P+S+1 > c s X p ) = 1 - / B- P"(|, %)dG( Vl , W ) ^ 0. □ 

Therefore, {X„}„>i satisfies condition P/(u„, a; n ), for any sequence of real vectors {u„}„>i and 
for any sequence {l n }n>i such that l n — > oo, corresponding to the multivariate version of Leadbetter's 
D-condition of local dependence (see, e.g., Leadbetter et al. Q~5] 1983). 

Now we shall see that {X„}„>i also satisfies the multivariate version of D" condition of Leadbetter 
and Nandagopalan ([IS], 1989). For a given sequence of real vectors {u n }„>i, we say that condition 
D"(u n ) holds if D(u„,ai n ) also holds and 

^E|=2 fen] Pfrl £ "n, X, < U„, X i+1 £ U n ) -> 

for some sequence {k n } n >i such that, as n — > oo, 

fc„ — s- oo, -> and k n ai n -t 0. 

Proposition 3.2 {X n }„>x satisfies condition D" (u„ ') . 
Proof Observe that 

P(Xl £ u£"\x, < u£ T) ,X i+1 £ u£ T) ) 
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Since each marginal sequence {X n _j} n >i satisfies condition D"(u^ 3 J) (Canto e Castro, [3] 1992), the first 
term above has null limit, as n — > oo. The second term above is upper bounded by, successively, 

£!<.,✓<* »E£a" ] p ( x hs > ufa\x iiS , < «<£',> < c s ,X i!S , V Y w ) 
= Ei<s, s '< d nT,^2 n] P(Xi, s > «fcUy < u<# < Y i+1 , s ,) 

= £i<,.,<^fe](i-f.(^ ) ))(i- /; , ff J 

< £ Ei< s ,^<^(l - F s (vi T Mn(l - Fy(«^)) - n(l - Mu£;',>/<v))) F : 

\ / F 8 t{u 8 , fC s /) 

which also converges to zero for any sequence k n — > oo, since by © we have n(l—F s i(u^)/c a i)) — > t*, > 0, 
as n — > oo. □ 



4 The multivariate extremal index and the domain of attraction 

of {X n } n >i 

A phenomenon also noticed in real data is that extreme events often tend to occur in clusters. The mea- 
sure that is used to capture the clustered extremal dependence is the extremal index (Leadbetter et al. 
|15| 1983). More precisely, the extremal index can be interpreted as the reciprocal of the limiting mean 
cluster size. A unit extremal index means no serial clustering and is a form of asymptotic independence 
of extremes. Figure [T] illustrates both clustering (symbol "•") and asymptotic independence (symbol "*") 
at high levels. 

The results of the previous section will allow us to compute the multivariate extremal index of {X„}„>i 
(Nandagopalan [3T] 1990). More precisely, if for all r € W\_ there exists normalized levels {u^ = 
(u^V, u n d d)}n>i, such that the sequence {rtP(Xi ^ u« )}„>i is convergent and D"(un ) holds, 
then {X„}„>i has multivariate extremal index if and only if, Vt 6 R^, sequence {nP(X 1 <ui T) ,X 2 £ 
ul 1 " ^)}n>i converges too. In this case, 

f(n,— ,r d ) — irnin-^oo (Xi ^ u ( T)) , t€K + , (10) 

(Ferreira, [7] 1994). In the sequel we shall use the copula function with notation 

C F ( Ul ,..,u d ) = F(Fi l (ui),...,F^\u d )), (ui, ...,Ud) 6 [0,l] d . 
Proposition 4.1 If F £ 13(H) then {X n }„>i has multivariate extremal index 

fl , , 1 logC gl (e-^,jeJ) 

f(n,-,r d ) = 1 - - — , — — t-, 

logC^e Tl ,...,e Td ) 

where I is the set of indexes in D for which Hj(x) = $> aj (x) = e~' T 3 , x > 0. Moreover, 

f 1 

J \ 1 r ; ' . // //, -!-„ 

is i/ie extremal index of {X n j} n >x, j = 1, 



G 



* * 



Figure 1: Marginal sample paths of a bivariate ARM AX process with unit Frechet innovations: the 
symbol "*" corresponds to c— 0.1 where the high values tend to occur almost singly and symbol "•" 
corresponds to c=0.8 with visible high values clustering. 



Proof By hypothesis, F € r D{H) 1 i.e., Fj £ V(Hj,{a n j > 0},{6 nj }), with Hj of the extremal type 



A, ^> a or <& Q and Cp(uV n , ..juV") — > Ch (ui, (ui, Ud) € [0,l] d . Thus we guarantee the 
existence of normalized levels Un = (u£ \ , u^ d J), for which condition D"(u^ r ' 1 ) holds. 
Moreover, for u^y = a n ,jHj 1 {e~ T: >) + b n j, j = 1, d, we have 

nP(X x £ u^) = n(l - FCan.iHf ^e-^) + 6„,i, .... a M f^ 1 (e _T *) + &m)) 

On the other hand, 

7iP(Xx < ui" r) ,X 2 £ ul T) ) = nP(X x < ui' r) ,cX 1 V Y 2 ^ u£ t) ) 
= nP(X x < ui T) )P(Y 2 £ uk x) ) = P(Xj < ulrV(l - G(u£i\ 



P(Xi<u^') 
P(X 1 <uk'°/c) 



■ (n(l - P(^, - n(l - P(^, 4^))) 



-> -logi7(i/f 1 ( e - T 0,...,^ 1 (e- rd ))+logH( J ffr 1 (e-^) I ---,^" 1 (e^)) 

= -\ogH(H^\e-^),...,H2 1 (er^)) + \ gH I {H^ 1 (er^),... 1 H~ 1 {e-<)) Il 

where I is the set of indexes in D for which r* given in (J5]) are positive, i.e., for which Hj(x) = $ Qj (x) = 
e _a: 3 , x > 0, and if/ denotes the marginal distribution of H corresponding to those indexes. Therefore, 
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applying (fTUf) . we have 

log^(ffr 1 (c-^),...,^- 1 (e--')) 

1 _ logC gl (e- T ',-,e-^*) 7 
logCjj(e- Tl , ...,e- r <i) 

Observe that if 7 = then 6(ti, r<j) = 1, Vt, and if 7 ^ 0, we have 

f l , ifjg£>-7 

leading to the assertion (|11[) . which corresponds to the univariate marginal extremal index already derived 
in Alpuim (pQ, 1989). □ 

The expression obtained for the multivariate extremal index function has the advantage of, once 
known/estimated the constants Cj and the marginal domains of attraction, we are only dependent on the 
attractor copula of {X n }„>i corresponding to the i.i.d. sequence with the same distribution F. 

Since we have Cp(uV n , uV") — > Ch(ui, Ud) uniformly in [0, l] d and Ch is continuous, we can 

n— »oc 

replace the discrete variable n by a continuous variable t and equivalently state t(l — C F (1 — x\/t, . . . , 1 — 
Xd/t)) — > log Ch (e~ Xl , e _a: ' i ), x € [0, oo) d . If we rewrite the result of the previous proposition as 



t— >oo 



(n, ...,T rf ) = 1 - lim 



f(l - C F/ (1 - TIC? 1 /*, . . . , 1 - TdC^/t)!) 



1 — lim 



t{X-C F [\-Ti./t,...,\-r d /t)) 

tp(u 3eI {F XlJ (x 1 ^>i-T 1 cr/t}) 



tP(UL l {Fx 1 AXij)>l-r 1 /t}) 



we can then estimate the multivariate extremal index through tail dependence functions estimators con- 
cerning Fj and F. For this issue see, e.g., Schmidt and Stadtmiiller ([22, 2006), Einmahl et al. ([5], 
2012) and references therein. 

Example 1 Consider F with F lt F 2 g £> '(A) andFj g X>($i), j = 3, d. IfC F (u u u d ) = exp(-(X)* =1 (- logu,) 7 ) 
7 > 1, i/ien -F £ T>(H). with Ch = C F , H± = H2 = A and Hj = $1, j = 3, d. Therefore, we have 

at \ 1 Ej=3( r J c i) 7 ) 1 ^ 7 / v _ w d 

0{r 1 ,...,T d ) = 1 ^ ^T77— > (n,-,r d ) e M +! 



(E- =1 r7)^ 



1 and 0j ■ = 1 — Cj, j = 3, d. 



Example 2 Consider F with F 3 e 2?($i), j = and C F (ui, u d ) = Aj=i u r Then F € 2?(77), 

TOtft, Cij = Cf, 77j = $1, j = 1, d. Therefore, we have 

\ 1^ 

6{<n, -,Td) = 1 - r ' r '- (ri,...,r d ) G R*. 
Vi=i 

The next result relates the domain of attraction of F with the one of G. 
Proposition 4.2 If F g V(H) then G g V(V) with V 3 - = H® 3 and 6j given in (17]). j g D, and 

v-y / l/#l l/8d\ 

(12) 
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Proof By hypothesis, Fj G V(Hj,{a n ,• > 0}, {&„,)}), j G D, i.e., F- l (a n .jXj + b n ,) — > HAxA and 

J ' ' n— foe 



C£(% /n , ..,uJ n ) C H {u 1 ,...,u d ), (wi,...,Wd) € [0,l] d . In addition, 

n— ►oo 



1 " fe) -if", '!>■.• 



From the stationarity relation in ([6]) , we have 



G~ (ci n jXj ~\~ b n j) y 



Therefore, 

//,;•<•,; • if//, (a,* q ! 
^"'Mzj) . if //.. -I- , . 

and thus Gj G T>(H^), j £ D, with 0j given in ([TT]). 
Now we look at the copula of G. We have 

F"(o nj ia;i + b n i,...,a n AXd + b nd ) -> #(xi, x d ) = C ff (Hi (xi ),..., H d (x d )) 

and 

* s '-' t„ - L F M s >-A*d 



= //'(-../ 

Again from the relation between F and G in we obtain 

G {a n ,iXi + b n ,i,...,a n ,dXd + On,d) -> — s- — • 

™ C^(i// :.r ; ;., /) 
Thus we can say that G G T>(V), where Vj = H^ 3 and 

...,x d ) = HA— ,j G /)^(o;i,...,Xd), 

or equivalently, 

C H (ui, ...,u d ) = C H (u\~ e \ ...,u d ~ d )Cv{u^,...,u e d d ), 
given 9j, j £ D, stated in (fTTl) . □ 

Observe that if 6j = 9, j = 1, . . . , d or if C# is the product copula then Cy = Ch- However, in general, 
relation (TT21 adds one more method to the existing ones of copulas construction (Joe 1997, Liebschen 
[30] 2008). Any MEV copula can appear in the limiting distribution H generated from this model, since 
H G V(H). Thus considering a MEV copula C and constants 0j g]0, 1], j = 1, . . . , d, we can derive new 
MEV copulas by applying the ratio rule in (fT2"j) one or more times. We shall illustrate the procedure by 
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considering that C is a Gumbel copula. More precisely, if C(m, Ud) = cx p(— l°g "j) 7 ) 1 ^ 7 ); 
7 > 1, then 

exp(-(£ti(-£l°g^) 7 ) 1/7 ) 



ex P (-(E J ti(-(i- 1 ) 1 °g^) 7 ) 1/7 ) 



is a MEV copula. The extremal coefficients, and thus the tail behavior, of this new copula present 
a greater variability of values when compared with the respective ones of the Gumbel. Indeed, for 
any J C D, we have Cj(u,...,u) = u € J with t C j = IJ] 1 ^. and C* r (u, ...,u) = u e J with ej" = 

5 Coefficients of tail dependence and tail independence 

Loosely speaking, tail dependence describes the limiting proportion of exceedances by one margin of 
a certain high threshold given that the other margin has already exceeded that threshold. The most 
used definition of tail dependence, provided in the monograph of Joe ([M], 1997), is the tail dependence 
coefficient (TDC): 

X = lim P(F Y (Y) > 1 - t\F x (X) > 1 -t). (13) 

We say that the random pair (X, Y) is tail dependent whenever A > and tail independent if A = 0. 

In the tail independent case, Ledford and Tawn ([T5JQ15] 1996/1997) proposed to model the null limit 
in (|13p by introducing a coefficient (77) to rule the decay rate of the joint bivariate survival function: 

P(F Y (Y) > 1 - t\F x (X) > 1 - t) ~ L{t)t 1/n -\ as 1 1 0, 

where L is a slowly varying function at 0, i.e. L(tx) / L(t) — > 1 as t 4- 0, for any fixed x > and 77 G (0, 1] is 
a constant. Coefficient 77 measures the degree of tail independence between r.v.'s X and Y. Observe that 
tail dependence occurs if 77 = 1 and L(t) -fr 0, as 1 1 0, and tail independence otherwise. The r.v.'s X and 
Y are called positively associated when 1/2 < 77 < 1, nearly independent when 77 = 1/2 and negatively 
associated when < 77 < 1/2. 

Both concepts can be naturally extended to a lag-r (r € No) formulation of a stationary d-dimensional 
sequence, {X„ = (X n j, X rly d)}n>i- More precisely, the lag-r TDC as 

A$(X) = ]imP(F r (X 1+r . r ) > 1 - t\F 3 {X^) > 1 - t) 
and the lag-r (r € No) Ledford and Tawn coefficient 77^, (X) defined by 

P(F r (X 1+rjl ) > 1 - t\Fj(X hj ) > 1 - i) ~ i(t)rj 1/r '^' (XK1 , as f 1 0, 

where i is a slowly varying function at 0. We denote Xjji(X.) = A^y(X) as the TDC between the jth 

and the j'th components, A^(X) = X^V (X) is the lag-r TDC within the jth sequence and A^)(X) is 
the lag-r cross-sectional TDC between the jth and the j'th sequences. An analogous description holds 
for the Ledford and Tawn coefficients, respectively, 77^/ (X), 77^ (X) and 77^, (X). 

In the sequel we shall denote wyt = 1 ( 1 — t). 
Proposition 5.1 //{X„} rl >i has stationary distribution then 

# <*>- 2 - us K 1 - (1 - * v***^ 

where Cjf denotes the common copula of (X n j , X n ji) , n > 1. Moreover, 



< 2 - 1 1 1 - C„, (1 - t, F, (c/^ t ))^fe J <2-i 1- jr^-T ) • ( 14 ) 
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Proof We have that 



hm ^&fl > 1 " h p y ( x ^f) > 1 - f ) 

Ho P(Fj{X hj )>l-t) 

2 - Urn i (1 - P(Fj (X hj ) <l-t, Fy (X w )<!-*)) 



2 - hm -Jl- P(F 3 (X hj ) < 1 - t,F f (X hj ,) < Fy(cJ, r wy t ))f[Gy(wy t / C y) 



Now the first result follows from ([5]). The second assertion is a consequence of the Frechet-Hoeffding 
copula bounds, i.e., max(ui + u<z — 1, 0) < C(u\, 1*2) < min(ui, 1/2), for all (ui, u-f) S [0, l] 2 . □ 



In the following we will state some consequences of this result, considering different situations for the 



domains of attraction of Fy . 



Proposition 5.2 Under the conditions of Provosition HOI we have 



^<X)=2-toI(l-C jj ,(l-M-4 <, ") I ^4 



if Fy e V(<P ajl ). Moreover, < A^)(X) < c-f '. 

Proof The result is straightforward since we have Fj>(cJ, r Wj>t) ~ 1 — tc™ J ' (see, for instance, Proposition 



3.3 in Ferreira and Canto e Castro OH 2008). □ 



whenever Fy has positive finite right end-point, with r£N. 



Proposition 5.3 Under the conditions of Provosition COl we have AvJ,(X) — and 77^/ (X) = 1/2, 



Proof Observe that FyiF^ll - t)c~ r ) ~ 1, and thus 

3 3 iin 



/|0 

i(l - C^(l - t, 1)(1 - t)) ~2 - |(1 - (1 - t) a ) ~*- □ 



Proposition 5.4 Under the conditions of Proposition 15. Jl we have A$(X) = and 1/2 < ^(X) < 

c -rk 

max(l/2,cj/ ) whenever Fy (cJ, r Wj>t) ~ 1 — £ 3 *' , /or fc > and r£N. 

Proof Just observe that the left and right hand-side of (fH| approximates, respectively, i and i + t c j' , 
as 1 1 0. □ 

Examples of d.f.'s satisfying Fy{cf, T wyt) ^ 1 — t 3 ' , fc > 0, include, e.g., Weibull of minimums (with 
d.f. F(x) = 1 — exp(— a^)) and Exponential (fc = 1). 

In a max-autoregressive context, the non negative associated tail independence (1/2 < 77 < 1) can 
also be described through a pARMAX process (Ferreira and Canto e Castro [TUIITT], 2008/10), whose 
logarithm corresponds to ARMAX. 

An illustration of the tail dependence between the marginals of {X„}„>i , for the three cases of domains 
of attraction can be seen in Figure O Observe that the copula's dependence is determinant: dependence 
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Figure 2: Scatter-plots for model (Xi, X2, X$) based on Gumbel's copula and marginals ARMAX, re- 
spectively, c = 0.8,0.1,0.1 and innovations distributed as unit Frechet, Exponential and Uniform; Left to 
right: points of (X\,X2), (Xx,Xs) and (X2,Xs) with Gumbel's copula dependence parameter 7 = 0.1 
(strong dependence) on the top and Gumbel's copula dependence parameter 7 = 0.9 (weak dependence) 
on the bottom. 

is evident whenever a strong dependent copula is used, whilst a weak dependent copula leads to an almost 
random scatter-plot. Figure [3] illustrates cross-sectional tail dependence of {X„}„>i, considering again 
the three domains. Observe the presence of some tail dependence for random pairs (Xj, Xj, ) whenever 
the lag-r apart j' th marginal is Frechet (first column plots) corroborating Proposition 15.21 An almost 
randomness can be seen in the other scatter-plots which is consistent with Propositions 15.31 and 15.41 



6 Marginal parameters estimation 

In this section we shall focus on the marginal ARMAX autoregressive parameter and the marginal tail 
index. 

The following result allow us to state an estimator for the ARMAX parameter Cj , j € D. 
Proposition 6.1 If Fqj and Gj are unit Frechet d.f 's, j £ D, then 

Cj = 2 -E^T) (15) 

where Xj is a r.v. with the stationary distribution of {X n j} n >i. 

Proof We shall use the result of Proposition 3.1 in Ferreira and Ferreira ([5], 2012). More precisely, 
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Figure 3: Cross-sectional scatter-plots for model (Xi, X 2 , X 3 ) based on Gumbel's copula and marginals 
ARMAX(0.5) with innovations distributed as, respectively, unit Frechet, Exponential and Uniform; 
Left to right and top to bottom: points of (Xi,x[ 2) ), (X U X^ 2) ), (Xi,x| 2) ), (X 2 ,x[ 2) ), (X 2l X (2) ), 
(X 2 ,xf } ), (X 3 ,X[ 2 } ), {X 3 ,X {2) ) e (X 3 ,X^ 2) ), where xj r) denotes the j th marginal lag-r apart. 
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considering F(x) = e~ x and s € N, we have 

E{F{X n ,Y) = E(F(c]X 0tj V VIU c^Y itj )>) 



s 



log^, yiii> ... >rnj) (c7",cJ" +i ,...,cjM) (16) 



J ' J ' " * " ' 3 

Ml „n— j 



Assuming that {A nj }„>i is stationary and taking limits in both of the members with s = 1, we have 



E(e~ x ^) = 1 ~ c 



leading to the assertion. □ 

As a consequence of this result, we verify that if {X n j} n >i is stationary than E(e~ x j ) € (1/2, 1). 
From (|15p we derive the estimator 

where = — EILi e _J<:i >J . Observe that no definite result can be obtained for Uj < 1/2, which may be 
an indication of an unsuitable model's choice. 

Based on ([3]), we have that X n j = V^i c *jY n -ij is the unique stationary solution of recursion (JTJ) 
(Davis and Resnick [5], 1989; Proposition 2.2). Therefore, an ARMAX process is ergodic (Stout [H], 1974; 

Theorem 3.5.8) and, since E(\e~ x j |) < oo, we have Uj —> E(e i ) almost surely (see, e.g. ergodic 
theory in Billingsley [5] 1995). Thus estimator cj is strongly consistent. The asymptotic normality is 
stated in the next result. 

Proposition 6.2 Under the conditions of Provosition [K7[ we have ^/n{Uj — e~ x i ) — > N(0,a 2 ) where 



a 



1 C J_ 



3 -2c, \2-c 3 J ^V(2-c i )(2-c J -cJ-4+ 1 ) K2 ~ c i 
Moreover y/n(6j - Cj) -> N(0, cr 2 (3 - 2c,)) . 

Proof The asymptotic normality also holds given the strong-mixing dependence structure with variance 
given by (Billingsley, 0, 1995) 

OO 

o- 2 = var(e~ x i 1 ) + cov((e~ x ^) (e~ x i+-)) . 

r=l 

Since Gj and Fqj are unit Frechet, according to the stationarity relation in ©, then 

3,-1 1 

= e I_c i , x > 0, j e D. 
For each reN, the joint d.f. of (X n j,X n+r j) is given by 

F ( x n j,x n+t , j )(x, y) = P(X nij < x A ycj r )P(Y n+1 tj < ycj r+1 , Y n+rJ < y Cj ), x > 0, j G D. 
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whose joint density, for y > xcj r > 0, is derived as 

9 9 „ , . d d ( -x- 1 —!— -fj -i r -i 



i=l 



— x — - — x u j —y - 



e — e , j G D 

x 2 (l-c 3 ) y 2 (l -Cj) 

Therefore, and after some calculations, we obtain 

poo rye ■ q £i 

E (e- x 7\- x 7 + 1 r)= I I 3 , ' , » --/' v. .v . (x,y)dxdy 



o Jo 



1 - c r 

j_ 



(2-c j ){2-c j - < f j -^+ 1 Y 

Now just observe that, by (fTTil) , 

E((e-V) 2 )= ^ 



2 



Considering g(x) — 2 — we have </^i?(e j )J = E[e i ) and, by the Delta Method, the 
second assertion holds. □ 

Other estimators can be found in literature. A strongly consistent estimator was earlier proposed in 
Davis and Resnick (0, 1989): 



= A 



x, 



i=2 1 1 



Another estimator, with a quite similar expression to our proposal, was derived in Lebedev (|17j. 2008). 

3,-1 i 

More precisely, for unit Frechet marginals Fj and Fqj(x) = Gj(x) = e , then 

Cj =2--, 
Pi 

with pj = P(X 2 ,j < X ld ) G (1/2, 1), and thus 

1 

Pi 



where pj = (n — 1) 1 Y^i=i <x t ■}• Note that a similar restriction to our method must be con- 

sidered, i.e., 1/2 < pi < 1. The consistency and asymptotic normality of this estimator can be seen in 
Ferreira (0, 2012). 

In what concerns the tail index ctj of each marginal j G D, it can be estimated through tail index 
estimators already stated in literature as Hill (in case otj > 0), Pickand's, maximum likelihood, moments 
or generalized moments estimator, whose asymptotic properties of consistency and normality still hold 
under an ARMAX dependence structure (Ferreira and Canto e Castro [10] 2008, Proposition 3.4). 
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